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ABSTRACT 

We study the SZ-effect-induced non-Gaussianity in the cosmic microwave background 
(CMB) fluctuation maps. If a CMB map is contaminated by the SZ effect of galaxies or galaxy 
clusters, the CMB maps should have similar non-Gaussian features as the galaxy and cluster 
fields. Using the WMAP data and 2MASS galaxy catalog we show that the non-Gaussianity 
of the 2MASS galaxies is imprinted on WMAP maps. The signature of non-Gaussianity can 
be seen with the 4 th order cross correlation between the wavelet variables of the WMAP 
maps and 2MASS clusters. The intensity of the 4 th order non-Gaussian features is found to 
be consistent with the contamination of the SZ effect of 2MASS galaxies. We also show that 
this non-Gaussianity can not be seen by the high order auto-correlation of the WMAP. This 
is because the SZ signals in the auto-correlations of the WMAP data generally is weaker than 
the WMAP-2MASS cross correlations by a factor / 2 , which is the ratio between the powers 
of SZ effect map and the CMB fluctuations on the scale considered. Therefore, the ratio of 
high order auto-correlations of CMB maps to cross-correlations of the CMB maps and galaxy 
field would be effective to constrain the powers of SZ effect on various scales. 

Key words: cosmology: theory - large-scale structure of the universe 



1 INTRODUCTION 

The thermal Sunyaev-Zel'dovich (SZ) effect, or the inverse- 
Compton scattering of the cosmic microwave background (CMB) 
photons by hot electrons, shifts the spectrum of the CMB photons 
to higher energy when the photons pass through the regions of cos- 
mic hot gas. It yields a CMB temperature T change at frequency v 
given by 



Ar sz (n) 
T 



( \( eX + 1 



where the dimensionless Compton y-parameter is 



y(n) = <tt dl 



,k(Te - T) 



(1) 



(2) 



where the integral is along the line of sight of n, and n = (/, b), 
b being the Galactic latitude and / the Galactic longitude, or is 
the cross section of the Thomson scattering and x = hv/kT. n e 
and T e are, respectively, the number density and temperature of hot 
electrons. Since groups and clusters of galaxies are hosts of hot gas 
and their distributions are non-Gaussian, the SZ effect will imprint 
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the non-Gaussianity of groups and clusters on the maps of the CMB 
temperature fluctuations (Cole & Kaiser 1988). 

Since the WMAP data became available (Bennett et al 2003a, 
2003b), many groups have detected the SZ effect signals with the 
cross-correlation between WMAP data and galaxy samples. The X- 
ray based catalogues of clusters and galaxies, including the North- 
ern ROSAT All Sky Galaxy Cluster Survey (Bohringer et al. 2000) 
and the ROSAT Brightest Cluster Sample (Ebeling et al. 1998) are 
found to be cross-correlated with the WMAP maps at the 2-5 a 
level (Hernandez-Monteagudo and Rubino-Martin 2004). The sam- 
ples of groups and clusters identified from the APM galaxy sur- 
vey (Maddox et al. 1990), and 2MASS (Jarrett et al. 2000) show 
cross correlation with W band data of the WMAP (Myers et al. 

2004) .The WMAP-2MASS correlation is also found by compar- 
ing it with model expectation (Afshordi et al. 2004). Using a semi- 
analytic model of the Intra-Cluster Medium, SZ signal was detected 
from 116 low redshift X-ray clusters at ~ 8a level (Afshordi et al. 

2005) . 

On the other hand, the WMAP data are found to be Gaussian, 
especially no significant non-Gaussian signals have been detected 
on scales of clusters (e.g. Komatsu et al. 2003). Therefore, a prob- 
lem is whether the positive results of detecting the SZ effect with 
WMAP-galaxy cross-correlation is consistent with the negative re- 
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suits of detecting non-Gaussian features of the WMAP data alone? 
In this paper, we try to reconcile the two results. That is, we want 
to show that 1.) the WMAP maps contain SZ-effect-induced non- 
Gaussianity; 2.) this non-Gaussian signal cannot be seen with the 
WMAP maps alone. 

We will study the non-Gaussianity induced by the SZ effect 
of the galaxy sample listed in the 2MASS extended source catalog. 
First, the SZ effects of 2MASS galaxies have been detected at about 
3 a level with the cross-correlation between the maps of WMAP 
and 2MASS (Myers et al. 2004; Afshordi et al. 2004). Second, 
the non-Gaussian features of the 2MASS samples have been exten- 
sively analyzed (Guo et al. 2004). It would be helpful to search for 
the WMAP's non-Gaussianity induced by the SZ effect of 2MASS 
galaxies. 

The paper is organized as follows. In §2 we present the data 
used for the cross correlation analysis. §3 analyzes the SZ effects of 
the DWT clusters, which are identified from the 2MASS galaxies 
with the discrete wavelet transfer (DWT) decomposition. The non- 
Gaussianity of WMAP induced by the SZ effect of the 2MASS 
DWT clusters is analyzed in §4. The discussion and conclusion are 
given in §5. Some math stuffs with the DWT algorithm are given in 
Appendix. 



2 SAMPLES AND ITS DWT DESCRIPTION 
2.1 Data of CMB temperature fluctuations 

We use the foreground cleaned WMAP maps, AT(n), of W and 
Q bands (Bennett et al. 2003a). The contamination of the galac- 
tic foreground is reduced by mask Kp2 (Bennett et al. 2003b). 
These maps were used for producing the WMAP first-year power 
spectrum of CMB temperature fluctuations. The foreground signal, 
consisting of synchrotron, free-free, and dust emission have been 
removed. 

The frequency of Q band is 40.7 GHz, and the SZ effect is 
(AT 3Z /T)q = -lMy, while the frequency of W band is 93.6 
GHz, we have (AT SZ /T)w = —1.56y. That is, at these two fre- 
quencies, the SZ signals AT SZ is less than 0. Therefore, the SZ sig- 
nals would be survived during the removal of foreground emission 
sources. 



2.2 2MASS-XSC galaxies 

We use the 2MASS extended source catalog (XSC, Jarrett et 
al. 2000), which covers almost the entire sky at wavelength be- 
tween 1 and 2 fim. The condition of selecting galaxies is taken 
to be K_m_k20fe, which measures the magnitude inside a ellipti- 
cal isophote with surface brightness of 20 mag arcsec~ 2 in K s - 
band. There are approximately 1.6 million extended objects with 
K a < 14.3. Most of the XSC sources at \b\ > 20° are galaxies 
(> 98%). The contamination mainly is from stars. The reliability 
of separating stars from extended sources is 95% at |fc| > 10°, 
but drops rapidly to < 65% at |6| > 5°. To avoid the contami- 
nant of stars, we use a latitude cut of |fc| > 10°. We also removed 
a small number of bright (K s < 9) sources by the parameters of 
the XSC confusion flag (ccjlag) and visual verification score for 
source (vc). They are identified as non-extended sources including 
artifacts. Moreover, to eliminate duplicate sources and have a uni- 



form sample, we use the following parameters: use_src = 1 and 
dup_src = . 

To select the range of K s , we use the standard log TV — log S 
test to examine the completeness of the sample. The number counts 
can be approximated by a power-law (Afshordi et al. 2004) as 



dN 
dm 



oc 10 K ' 



(3) 



The XSC sources with |6| > 30° and 12 < K s < 13.7 are believed 
to be galaxies with 99% reliability (Mcintosh et al. 2003). For this 
sample, the index k is found to be 0.641±0.006. If considering this 
ft to be the standard, the completeness of a sample can be estimated 
by the deviation of (d7V/dm) samp i c from the standard, i.e. 



C(m) 



(dN/dm) 



(4) 



(dN/dm) sample 

where the standard sample is taken to 12 < K s < 13.7 and 
\b\ > 30°. It has been shown by Guo et al. (2004) that the com- 
pleteness C(m) is equal to 1 for sample of |6| > 10° in the range 
11 < K s < 13.7. The factor C(m) is obviously larger than 1 for 
sample of |6| > 30° when K s < 10.0. This indicates the cata- 
log to be contaminated towards the bright end. On the other hand, 
C(m) drops below 0.9 when K 3 > 14.0. Thus we use a cut of 
10.0 < K a < 14.0 to ensure our sample to be complete greater 
than 90%. This sample contains 987,125 galaxies with median red- 
shift z ~ 0.1. It gives a 2-D number density field of galaxies, 
Ps( n )- 

2.3 DWT variables of maps 

To analyze the cross correlation and non-Gaussian signals, we will 
use the variable given by the discrete wavelet analysis (DWT) de- 
composition. That is, statistical properties are measured by DWT 
mode-mode correlation, and non-Gaussianity is detected by high 
order statistics of the DWT variables. Since the DWT mode is lo- 
calized in phase (position and scale) space, and has regular shape, 
it can also be used for 1.) handle the bad pixels in maps; 2.) identify 
clusters from galaxy map. 

Since the SZ effect is on small scales, one can subject the 
CMB temperature maps to an equal-area projection by the Lam- 
bert azimuthal algorithm: 



(5) 



Ry/2 - 2 1 sinfc| cos /, 
X2 = Ry 2 — 2| sin b\ sin /, 

where R is a relative scale factor,/? is the Galactic latitude and / is 
the Galactic longitude. This hemisphere scheme projects the whole 
sky into two circular plane, northern and southern sky. 2MASS 
galaxies are also described by this format. We select a square with 
123°. 88 x 123°. 88 in the central part of each circular plane. We 
have two fields of 123°. 88 x 123°. 88 in northern and southern 
sky. Both the WMAP map and 2MASS galaxy distribution are fully 
overlapped with each other. We will use coordinate x = (xi,X2) 
to replace n = (I, b) below. The 2-D maps of AT(n) and p g (n) 
will be written as AT(x) and p s (x). 
The DWT variables are defined as 



AT: 



j.l 



/ 0j,i(x)dx 



AT(x)0 j]1 (x)dx, 



(6) 



The notations of the 2MASS parameters used in this para- 
graph are from the list shown in the 2MASS Web site 
http://www.ipac.caltech.edu/2mass/releases/allsky/doc 
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pj.i 



Pfl( x )<Ai,i( x ) dx , 



and 



r 3 
c j-l 



AT(x)^j,i(x)dx, 



Ps( x )V'j,i(x)dx, 



(7) 



where 0j,i(x) and ^>j,i(x) are, respectively, the scaling func- 
tion and wavelet. For our 2-D samples, the DWT variables 
A3j i and Pj,i describe, respectively, the mean temperature and 
the mean number density of galaxies in the cell (j,l), which 
has size 123°.88/2 J1 x 123°.88/2 J2 and at position around 
[/i(123 .88)/2 J1 ,; 2 (123 .88)/2 J1 ], where ji and j 2 can be any 
integral, and h = 0, ...2» _ \ h = 0, ...2 J ' 2_1 . Thus, the DWT 
index j corresponds to an angular scale of 123°. 88/2-'. The an- 
gular distance between modes 1 and 1' at scale j is given by 
6 = 123.88° 1 1 - 1' | /2 3 . The wavelet variables (WFCs) and 
describe, respectively, the fluctuations of temperature and galaxy 
density on scale j and position 1. Some details of the algorithm 
with the DWT is given in Appendix A. 

The 2MASS XSC galaxies are resolved to 10 . Our analysis 
of the 2MASS sample can reach to angular scale of about 0.01 de- 
gree. However, on scales less than 9 = 123°.88/2 9 = 0°.24, the 
WMAP data are dominated by noise. In calculating the WMAP- 
2MASS cross correlation or WMAP's auto-correlation, we will use 
only scales of j = 8. 

Since the DWT variable ATjj is localized, it can be used to 
handle the foreground masks. A standard technique to treat con- 
taminated pixels is zero-padding (Pando & Fang 1998). That is, 
1.) put zero data at the masked pixels, and 2.) off-count the DWT 
modes (j, 1) located at the masked pixels. We test the zero padding 
by simulation samples generated with code HEALPix 2 . The results 
show the Kp2 masked samples with zero padding yield the same 
statistical properties as the original simulated samples, at least, up 
to the 4 th order. 

We should emphasize that the Lambert projection will vio- 
late the rotational invariance on the spherical surface. We will not 
use the rotational invariance in our analysis. Moreover, we show in 
Appendix B that the Lambert projection does not cause false non- 
Gaussian features. That is, if the original map is free from non- 
Gaussian correlations, the projected map is also free from these 
correlations. The Lambert projection is legitimate for our DWT 
analysis. 



3 THE SZ EFFECT OF 2MASS CLUSTERS 
3.1 DWT clusters of 2MASS galaxies 

The SZ effect is sensitive to hot gas clouds. If the mass density 
and temperature of gas are proportional to the number density of 
galaxies, one can identify hot gas clouds on scale j by modes (j, 1) 
with the high These modes (j, 1) are called DWT clusters on 
scale j. The general method of identifying DWT clusters with pj.i 
have been studied with simulation and real samples (Xu et al. 1999, 
2000). It showed that the clusters identified by top pj,i statisti- 
cally are the same as the clusters identified by the friend-of-friend 
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Figure 1. Cross correlation (T(|Z — [eq.(8)] between the WMAP maps 
of W (top) and Q (bottom) bands and top 500 2MASS DWT clusters on 
scale j = 8. The angular scale of \l-l'\ is \l - «'|123.88/2 8 degree. 



method if the mean size of the friend-of-friend identified clusters is 
the same as that of DWT clusters. 

The scale of DWT variables is well defined, and therefore, it 
is parameter-free. On the other hand, the friend-of-friend algorithm 
needs the so-called link parameter. These parameters may introduce 
uncertainty in the correlation analysis. Moreover, The clusters iden- 
tified by the friend-of-friend method usually have very irregular 
shapes (Jing & Fang 1994), it is inconvenient to estimate the statis- 
tical significance of the cross correlation between CMB maps and 
galaxy catalog. On the other hand, the variables pjj and ATjj are 
in the same spatial cell of DWT mode, the statistical significance 
with pjj and ATj,; is unambiguous. The DWT scaling functions 
are orthogonal from each other, different DWT clusters consist of 
different galaxies. This is also useful for statistical analysis. 

We identified, in this paper, the DWT clusters on j = 8, which 
corresponds to angular scale 123°.88/2 8 ~ 0°.5, or length scale 
~ 1.8 h -1 Mpc at the median redshift of the sample, which is the 
scale of clusters of galaxies. We selected 500 cells with top pjj, 
which are the peaks with 6.4<r and higher, a being the variance of 
the fluctuations of the number density field of galaxies. The number 
of galaxies in the top 500 cells is from about 20 to 60. Although the 
500 cells is identified with 2-D sample, they most likely contain 
projected clusters or groups. 



3.2 The contamination of SZ effect 

We perform the cross-correlation between the WMAP data and the 
2MASS DWT clusters by 



AT(|/ — l'\) = (CjjAXjj/), 



(8) 



where the variable Cjj is taken to be 1 for mode (j, 1) correspond- 
ing to a DWT cluster on scale j = 8, and Cj,i = 0, other modes. 
The average (...) overs all possible \l — 1'\ of the sample. Therefore, 
AT(\l — l'\) is an average CMB temperature fluctuations with a 
distance \l — l'\ from DWT clusters on scale j — 8. 



4 Cao, Chu & Fang 



Figure 1 presents the cross correlation AT(\l — l'\) ofW and 
Q band maps with the top 500 2MASS DWT clusters. For each 
DWT cluster, one can have a correlation AT(\l — Z'|). The solid 
line and error bar in each panel of Fig. 1 are, respectively, given 
by the average and 1-a variance of AT(\l — l'\) of the consid- 
ered 500 top DWT clusters. Figure 1 shows anti-correlation of the 
DWT clusters with AT at \l — l'\ = with temperature decrease 
AT 3Z = AT(\l - (AT) ~ -15 ± 10 /iK, or the Comp- 

ton {/-parameter to be ~ (3.7 ± 2.4) x 10 -6 . This result is about 
the same as that given by 500 2MASS galaxy clusters selected by 
friends-of-friends algorithm (Myers et al. 2004). The frequency- 
dependence \(AT SZ /T) Q \ = 1.91i/ > \(AT SZ /T) W \ = 1.56j/ 
didn't shown in Figure 1. It is probably because the difference be- 
tween (AT az /T)q and (AT 3Z /T)w is much less than the errors 
of Figure 1. 

We also can see from Figure 1 that for Q band also have a 
weak anti-correlation at \l — l'\ = 1 with the level of AT SZ ~ 
— 10 ± 7 p K. This result indicates that the distribution of hot gas 
probably is not simply proportional to the number density of opti- 
cal and infrared galaxies, but is more spread than the distribution of 
galaxies. Recently, cosmological hydrodynamic simulation shows 
indeed that the hottest gas generally distributes more spread than 
galaxies (He et al. 2005). Therefore, it would be reasonable to con- 
sider that hot gas is spread in the range of \l — l'\ < 1 around a 
DWT clusters. 

We did not found significant anti-correlation with the DWT 
clusters more than top 500, or on scales larger than j — 8. There- 
fore, the 500 top DWT clusters give an estimation of the contam- 
ination on the WMAP maps induced by the SZ effect of 2MASS 
galaxies. On scale of j = 8, the map has divided into 2 x 256 2 
cells. The 500 j — 8 DWT clusters do not overlapped from each 
other, and therefore, there are about 500/(2 x 256 2 ) = 0.4% area 
of the temperature maps of the WMAP to be contaminated by the 
2MASS SZ effects with the order of AT SZ ~ -15 ± 10 p K. Con- 
sidering the Q band signal at \l — l'\ = 1, the contaminated area 
would be about 1%. 



3.3 Mock samples of SZ contaminated CMB maps 

If the anti-correlation shown in Fig. 1 is due to the SZ effect of 
2MASS clusters, we can mimic a SZ-effect-contaminated CMB 
map by the following mock sample 



AT(x) = AT cmb (x) + AT sz (x) + AT second (x), 



(9) 



where AT cm b(n) is the primeval temperature fluctuations, and 
AT secon( j(n) is due to secondary effects other than the SZ effect, 
such as the ISW effect and microwave point sources. The term 
AT sz (x) of eq.(9) is given by eqs.(l) and (2). We consider only 
the SZ effects caused by hot electron in 2MASS galaxies, which 
are in a small redshift bin around z ~ 0.1, it is described by sur- 
face density p 9 (x). Therefore, we rewrite the Compton parameter 
as y ~ (a T /m e c 2 )n c ° l (n)T e (n), where n c ° l {a) is the column 
density of electrons, and T e (n) is the density- weighted mean of 
temperature of electrons in the redshift range of 2MASS. Gener- 
ally T e ^> T cm b- Therefore, the thermal SZ effect can be estimated 
as AT sz (x) oc n= oi (x)T(x). 

If galaxies trace hot baryon gas, we have approximately 
nl°\x.) oc p 9 (x). The relation between temperature T and mass 
density p g actually is complicated, because the cosmic baryon gas 
is multiple phased (e.g. He et al. 2004). For a given dark matter 
(or baryon matter) mass density, the PDF (probability distribution 
function) of the temperature of baryon gas covers a large range 
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Figure 2. Cross correlation (T(\l — l'\)) [eq.(8)] between mock sample A 
and top 500 2MASS DWT clusters on scale j = 8. The angular scale of 
\l - /'| is \l - Z'|123.88/2 8 degree. 

from 10 4 to 10 6-7 K. Nevertheless, the relation between the mean 
temperature and mass density of the gas can be approximated as 
a polytropic relation T e oc p 9 ~ , where a is about 1.5 (He et al. 
2004). Thus, we have AT sz (x) oc p g (x). Therefore, to mimic the 
thermal SZ effect, it would be reasonable to assume 



AT sz (x) = -/ 



((^)2)l/2 



P 9 (x), 



(10) 



where ej x and e| " are, respectively, the wavelet variables (WFCs) 
of AT cm b(x) and p" (x) [eq.(7)]. Subjecting eq.(10) to a wavelet 
transform [eq.(7)], and considering AT >• AT SZ , AT SOCO nd, we 
have 



f=i(^) 2 ) l/2 /mr) 



-Ts2 N l/2 



(ID 



Because ((e^i) ) and ((??") ) are, respectively, the pow- 
ers of the fields of CMB temperature fluctuations and SZ effect on 
scale j (Fang & Feng 2000), the parameter / is the ratio between 
the two powers on scale j. 

Thus, with eq.(9) we can construct mock sample of SZ-effect- 
contaminated CMB maps by the following steps. First, the term 
AT cm b(x) in eq.(9) is produced by the HEALPix simulation. Sec- 
ond, the term AT sz (x) is given by eq.(10), in which p(x) is taken 
to be the value given by the 2MASS map if the position x is within 
the cells of the top 500 clusters, otherwise, p 9 (x) is taken to be 
equal to zero. We ignore the term AT socon d(x). Finally, we esti- 
mate / by doing cross correlation between the mock sample eq.(9) 
and 2MASS DWT clusters. Figure 2 plots AT(\l-l'\) vs. \l-l'\ for 
/ = 0.01 0.1 and 0.2. The best fitting to the observed SZ tempera- 
ture change AT 3Z ~ — 15±10 /iKis given by / ~ 0.1, which is ra- 
tio of the powers of SZ effect to the CMB fluctuations. This result is 
consistent with the estimation given by semi-analytical model and 
simulation (e.g. Cooray et al. 2004). This sample is called Mock A. 

Unlike Figure 1 , Figure 2 does not show the tail of AT ( j / — / ' | ) 
at \l — l'\ ~ 1. This is because we assume n e (x) oc p 9 (x), i.e. no 
hot gas locates out side of a DWT cluster. If we assume that hot 
electron n e (x) can exist not only at the cells of DWT clusters, but 
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Figure 3. Cross correlation (T(\l — l'\)) [eq.(8)] between mock sample B 
and top 500 2MASS DWT clusters on scale j = 8. The angular scale of 
\l - l'\ is \l - /'|123.88/2 8 degree. 



Figure 4. C. ' (|1 - l'|), vs |1-1'| atj = = (8,8) for CMB simu- 

lation samples. The error bars are the range of 90% of 100 samples. 



also in their nearby cells, we have sample of Mock B. The cross 
correlation between mock sample B and 2MASS DWT clusters are 
shown in Fig.3, which does show a tail of AT(\l — 1'\). We see that 
the sample of / = 0.1 is also basically the same as observation. 

The term AT SOCO nd (x) of eq.(9) comes from the ISW and mi- 
crowave point sources. The ISW effect is mainly from potential 
with linear evolution. The ISW effect given by nonlinear evolution 
is very small (Seljak 1996; Tuluie et al. 1996). One can ignore this 
effect if we focus on non-Gaussian behavior. If a microwave point 
source is from the 2MASS galaxies, their contributions to the SZ ef- 
fect have already been included in the WMAP-2MASS cross corre- 
lation. If the microwave sources are not from the 2MASS galaxies, 
they generally are uncorrelated with 2MASS distribution p 9 (x), 
and therefore, the SZ signal of the 2MASS-WMAP cross correla- 
tion would not be hurt by this sources. 



4 NON-GAUSSIANITY INDUCED BY SZ EFFECT OF 
2MASS GALAXIES 

4.1 Non-Gaussianity detectors 

Using the mock samples, we try to estimate which non-Gaussian 
features are detectable with the WMAP samples. Effective detec- 
tors of non-Gaussianity are given by high order auto-correlations 
and cross-correlations of wavelet variables (WFCs) AT cm t,(x) and 
Pg(x) (e.g. Pando et al. 1998, 2002; Guo et al. 2004; McEwen et 
al. 2004; Cruz et al. 2005). The normalized high order correlation 
of wavelet variables is defined as 



C?' 2 (|l-l'|) = 



Cf«(|l-l'|) 



((^■i) p (gj,i<n 
<(Aj,i)*><(4i')«> ! 



(12) 



: 0, 



where Aj } \ and Byy can be ef t or e? Because (ef \) 
the number p and q of eq.(12) should be even integer. When A^i = 
Bj,i, C?' 2 (0) is the kurtosis of the field considered. Therefore for 
a Gaussian field, we have 



3 \l-l'\=0 
1 \l-l'\>0 



(13) 



In Figure 4, we plots Cj^fll - l'|) for CMB map AT cmb (x) 
produced by the HEALPix simulation. As expected, this sample is 
Gaussian. 

On the other hand, the sample of 2MASS-XSC galaxies 
is non-Gaussian (Guo et al. 2004). Figure 5 plots the \l — l'\- 
dependence of C?' 2 (|l - l'|) for the 2MASS galaxies. We see that 
Cj 2 ' 2 (0) is significantly larger than 3, while at all other points, 
i.e. |1-1'| / 0, we have C?' 2 (|l-l'|) = 1. That is, the non- 
Gaussianity of 2MASS samples measured by C?' 2 (|l — l'|) is lo- 
calized. It depends only on the density distribution of galaxies in 
the area considered, regardless galaxies in other places. With this 
feature, one can say that the mock sample eq.(10) contains all the 
2MASS non-Gaussian features if they are measured by detectors 
eq.(12). Of course, if consider the contribution of hot electron in 
nearby areas of DWT clusters, like in mock sample B, the detector 
C 2 ' 2 (|l- 1' | ) would also show a tail till to about \l-l'\ = 1. How- 
ever, the non-Gaussianity measured by C- ' (|1 — l'|) still is local- 
ized in the sense that the hot electron distribution is determined by 
the position of the DWT clusters considered, regardless galaxies in 
other area. 



4.2 Signature of non-Gaussianity induced by SZ effect 

To detect the WMAP non-Gaussianity given by SZ contamination 
of 2MASS galaxies, we use the 4 th order cross-correlation between 
the CMB maps and 2MASS galaxy distribution defined as 



/"2.2m ((d) 2 ( g j,i') 2 ) 

(|1 l|)c r oss- ((eTi)2)(( ^ )2r 

Using eqs.(9) and (10), we have 

Cf 2 (|l-l'|) 0rO ss = 



(14) 



(15) 
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Figure 5. C 2 ' 2 (|l 



3 l'|), vs |1-1'| atj = 

galaxy samples of north (solid) and south (dashed). 



for 2MASS 



(( 



~T r „ 



Hi 3 ,) 2 ) 



((efj) 2 ^,) 2 ) 



<(ei,) 2 ><(^) 2 > <(?r,?) 2 ><(^) 2 > 

Since there is no correlation between primeval CMB map and 

galaxies, we have <(ejr b ) 2 (e?,i>) 2 > = ((^f)^) 2 )- 
Therefore, the first term on the right hand side of eq.(15) is always 
~ 1. It is irrelevant to the non-Gaussian features. The non-Gaussian 
signal fully comes from the second term on the right hand side of 
eq.(15). 

If the galaxy field is Gaussian, we have 

1 + 3/ 2 
1 



Cf 2 (|l-l'|)c 



1.03 



l'\ = 
l'\ > 



(16) 



duce C 2 ' 2 (0) 



The cross-correlation Cj : 



where we used / = 0.1. Thus, we may expect that the non- 
Gaussian kurtosis of 2MASS galaxies shown in Figure 5 will in- 
cross > 1.03. 

1 2 "1 - l'|) cross for the WMAP map 
and 500 2MASS DWT clusters is shown in Fig. 6. In the top panel 
of Fig. 6, the solid and dotted line are for W and Q band respec- 
tively. The black points and error bars are the mean and 90% range 
of 100 simulation samples without SZ term. Therefore, Fig. 6 do 
show C J 2 ' 2 (0) C ross > 103. The middle and bottom panels of Fig. 6 
are, respectively, the Cj' 2 (|l — l'|) C ross vs. |1 — l'| for mock sam- 
ple A and B with / = 0.1. The error bars are the 90% range of 100 
mock samples. From Fig. 6, one can conclude that 1.) There are 
positive signals of the 4 th cross correlation between the wavelet 
variables of WMAP data and 2MASS clusters; 2.) The 4 th cross- 
correlation is consistent with the estimation of mock samples A and 
B, and therefore, this non-Gaussian feature probably is from the SZ 
effect of 2MASS galaxies. 

4.3 Cross-correlation vs. auto-correlation 

The last problem we should study is whether the SZ-effect- 
caused non-Gaussianity can also be seen with the high order auto- 
correlations of WMAP maps. This point can be analyzed with the 
4 th order detector defined as 



Figure 6. C?' 2 (|l - l'|) C ross between a.)(top) WMAP map and top 500 
2MASS DWT clusters on scale j = 8, in which solid and dotted line are for 
W and Q band respectively, and the black points and error bars are the mean 
and 90% range of 100 simulation samples without SZ term; b.)(middle) 
the mock sample A and top 500 2MASS DWT clusters on scale j = 8; 
c.)(bottom) the mock sample B and top 500 2MASS DWT clusters on scale 
j = 8. The error bars of b.) and c.) are the 90% range of 100 mock samples 



c 



|l-l'|)cmb = 



<(^,) 2 )<(^) 2 >' 



(17) 



This detector actually is similar to (7. 
instead of e? j by ef } . 

Using eqs.(9) and (10), we have 



2 (|l-l'l) 

cross [eq.(14)], but 



2 ' 2 (|l-l'|)cmb = 

terms irrelevant to nongaussianity 



(18) 



+/ 



4 <(g|,T) 2 (g|,?0 2 > 

<(e?J) 2 )<(e?J,) 2 >' 



Comparing eq.(18) with eq.(15), it is clear that the non-Gaussian 
term in eq.(15) is proportional to / 2 , while in eq.(18) to / 4 . There- 
fore, the SZ signal in C?' 2 (|l — l'|) C mb is weaker than that in 



C? ,2 (|l — l'|) cross by a factor / 2 , which can be as small as ~ 10~ 2 
if considering mock samples A and B. Thus, it is not surprised that 
this non-Gaussianity is undetectable with current CMB maps like 
WMAP. 

Figure 7 gives C 2 ' 2 (|l - l'|)cmb vs. |1 - l'| for W and Q 
bands of WMAP data, and for mock samples A and B. It shows no 
significant non-Gaussianity. Therefore, the detector of cross cor- 
relation C?' 2 (|l — I'D cross is effective to pick up the SZ-effect- 
induced non-Gaussianity, while the detector of auto-correlation 
C?' 2 (|l — l'Dcmb is n °t sensitive to that non-Gaussinaity. This re- 
sult is not limited with the 4 th order detectors, but it should be a 
common feature of high order detectors eq.(12). To replace (e? j) 2 
with (?j^i) 2 will lead to a factor f 2 in the non-Gaussian term, and 
therefore, the cross-correlation detectors generally are more sensi- 
tive to the SZ-effect induced non-Gaussianity than auto-correlation 
of CMB maps. 
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Figure 7. C? ,2 (|l - l'|) cmb vs. |1 - l'| for W (top) and Q (bottom) bands 
of WMAP data. The mock samples of A and B are also shown. The error 
bars are the 90% range of 100 simulation samples. 



5 CONCLUSION 

With the DWT clusters of 2MASS sample, we confirmed the pre- 
vious results of the SZ effect of 2MASS galaxies on the WMAP 
data (Myers et al. 2004; Afshordi et al. 2004). For the foreground 
cleaned WMAP maps of W and Q bands, there are about 0.4-1% 
area of the maps is contaminated by the SZ effect of 2MASS DWT 
clusters with the Compton parameter y — (3.7 ± 2.4) x 10 -6 . 
With mock samples of the 2MASS SZ effect, we show that the non- 
Gaussianity of the 2MASS galaxies is imprinted on WMAP maps. 
This non-Gaussianity can be seen with the 4 th order cross correla- 
tion between the wavelet variables of the maps of the WMAP data 
and 2MASS clusters. The intensity of the 4 th order non-Gaussian 
features is consistent with the estimation given by the mock sam- 
ples of the SZ effect of 2MASS galaxies. We show also that this 
non-Gaussianity can not be seen by the 4 th and higher order auto- 
correlation of the WMAP maps. 

The space-scale decomposition of the DWT is powerful tool 
for studying the SZ effect. The DWT variable can be applied to find 
clusters from galaxy samples, to identify the SZ effect of the DWT 
clusters, to detect the non-Gaussianity with high order correlation. 
With this method, all statistics are based on the same set of the 
DWT variables with well defined scale, position and shape of the 
modes, the relation between the SZ effect contamination and non- 
Gaussianity of the CMB maps can be measured without ambiguous 
parameters. 

We also found, because the ratio between powers of SZ maps 
to the CMB maps generally is much less than 1, the cross cor- 
relation between the WMAP map and galaxy distribution would 
be more effective to pick up non-Gaussianity than auto-correlation 
of the WMAP map. Due to the limitation of the resolution of 
WMAP data, we studied only the SZ effect on scale j = 8, 
or ~ 0°.5. Once higher resolution map is available, the ra- 
tio C?' 2 (|I - l'|) cm b/C?' 2 (|l - l'|)cros S would be able to provide 
more information of the ratio of the powers of SZ effect and CMB 
fluctuations on different scales. 
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APPENDIX A: ALGORITHM WITH THE DISCRETE 
WAVELET TRANSFORM (DWT) 

Consider a 1-D density fluctuation 5(x) on a spatial range from 
x — to L. We divide the space into 2 3 segments labelled by I = 
0, 1, ...2 3 — 1 each of size L/2 3 . The index j is a positive integer 
and gives the length scale L/2 3 . The larger the j is, the smaller the 
length scale. Any reference to a property as a function of scale j 
below must be interpreted as the property at length scale L/2 3 . The 
index I represents position and it corresponds to the spatial range 
IL/2 3 < x < (I + l)L/2 3 . Hence, the space L is decomposed into 
cells 

The discrete wavelet is constructed such that each cell (j, I) 
supports a compact function, the scaling function (j>j^(x). In our 
calculations, the Daubechies 4 wavelet (Daubechies, 1992) are 
used. The scaling function satisfies the orthogonal relation 



8 Cao, Chu & Fang 



1 



(Al) 



where S K is Kronecker delta function. The scaling function <f>j : i (x) 
is a window function on scale j centered at the segment /. The 
normalization of the scaling function is J 4>j tl (x)dx = (L/2 J ) 1/2 . 
For a field p(x), its mean in cell (j, I) can be estimated by 



C 3,i 



Jo 4>j,i( x )dx 



where ; is called scaling function coefficient (SFC), given by 



p(x)<j>j,i{ x )dx. 



A 1-D field p(x) can be decomposed into 

23-1 

P^)=^2e p jtl 4> i ,i(x) + 0{>j). 



(A3) 



(A4) 



The term 0(> j) in eq.(A3) contains only the fluctuations of the 
field p(x) on scales equal to and less than L/2 j . This term does not 
have any contribution to the window sampling on scale j. 

The fluctuation on scale j in cell (j, I) is given by the WFC 
defined as 



e j,i = I P{x)i>j,i( x ) dx - 

10 



The wavelet ipjj (x) are orthogonal with respect to j and i 



/ 



i>j : i{ x )i>i',i'i x )dx = 8^i,5f tj ,, 



They are also orthogonal with <j)jj (x) as 
<Pj',l'( x )tpj,l( x ) dx = if 3 < j- 



(A5) 



(A6) 



(A7) 



The orthonormality eq. (Al) insures that the set of Pjj or ; 
do not cause false correlations. When the "fair sample hypothesis" 
(Peebles 1980) holds, the average over the ensemble of the random 
field can be estimated by the average over modes (j, I). 



APPENDIX B: LAMBERT PROJECTION AND GAUSSIAN 
FIELD 

Bl White power spectrum 

For a field on spherical surface T(8, (f>), we have 

T(fi) = ^a im y ;m (fi), (Bl) 

im 

where Y !m (f2) is spherical harmonic function. The coefficients a; m 
are given by 



air, 



T(fi)Y, m (fi)dfi, 



(B2) 



where dQ = sin 6d6d<j>. For a Gaussian field T(fi), we have 
{aimdi'm') = <H5u'5 mm /, (B3) 

where (...) means the average over the ensemble of samples. 5u> 
means no scale-scale correlations. 



Consider a coordinate transfer, such as the Lambert projection 
eq.(3), as 



xi = Xi(fi) = Rsjl - 2| sin0| cos0, 
x 2 = X 2 (Q) = Ry/2 - 2 1 sin0| sin0, 
the field T in coordinate (xi , X2) is then given by 



(B4) 



(A2) F{ Xl ,X2) = J 8 D [x 1 -X 1 {Q)]8 D [x 2 -X 2 (n)]T(n)dn, (B5) 



where S D is the Dirac delta function. Thus, the wavelet coefficients 
of the field F(xi, X2) are 



/ 



F(xi,X2)lpjp(xi, X2)dX\dX2 



(B6) 



where ipj P (xi, X2) is 2-dimension wavelet basis on space (xi, £2). 
(In order to avoid the confusion between the index / for wavelet and 
spherical harmonic, we use p for the position index of wavelet in 
this Appendix). 

From eqs.(B5) and (B6), we have 

{e jp e j/p ,) = (^J i> jv {X 1 {Q),X2{Q.)) (B7) 

ip fpl (Xi(n'),X2(n'))X2(n'))T(Q.)T(n')dndQ') = 
J i, jp (x, (n) , x 2 (n))^- v (^i ( n ')> x 2 (n')) 

y^y^(aimai'm')yim(^)Y;>m'(^')rf^^ / - 

For a Gaussian field eq.(B3) with constant of, we have 

(e jp e j/p ,)oc J Vip(Xi(n),X 2 (n)) (B8) 

V J -v(Xi(n , ),x a (n , ))5^y, m (n)Y, m (n , )<«Mn'. 



Considering 

2 Y im (n)Y im (n') = <5 D (n - n'), 



(B9) 



we have then 



V> « (BIO) 

J i> 3P (x 1 (n),X2(n))^ jlp ,(x 1 (n),X2(n))dn, 

For the equal area projection of the Lambert transform eq.(3) or 
eq.(B4),we have 



dQ = RT 2 dXxdX2. 
Thus, eq.(BlO) gives 

ipj p (Xi, X2)ipj> p > (Xi, X2)dXidX2 — 8jji5 pp i, 



(BID 
(B12) 



The last step is based on the orthogonal-normal condition of 
wavelet. Therefore, {ejiijn 1 ) is diagonal in terms of scale index 
j, i.e. free from scale-scale correlation. 

B2 Non-white power spectrum 

If (a im a ; / m /) = of 8w8mm' and af is not constant, the covariance 
(ijpij'pi) will still be diagonal or quasi-diagonal with respect to 
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In this case, we should consider the locality of wavelets, i.e. 
in the Fourier space, the wavelet tpji is localized in the area around 
scale j. From eqs.(B5) and (B6) we have 



ipj P (x 1 ,x 2 ) 



(B13) 



8 D [ Xl - Xi(n)]5 D [:r2 - X 2 (0.)]T(Q.)dQdxidx2. 
Using eq.(B9) and (B2), eq.(B13) yields 



ipj p (xi,X2)S D [xi — Xi(Q)] 



I 

8 D [x 2 - X 2 (n)]S D (Q. - Q')T(Q')dQ.'dQdxidx2 
^]a !m / ipj P (xi,X2)S D [xi — Xi(Q.)] 

l,m •* 

8 D [x 2 - X 2 {tt)]Yi rn (Q)dVLdx l dx2 



(B14) 



where 



ip jp (lm) = J ij jp (x u X2)5 D [xi -Xi(fi)] (B15) 

S D [x 2 - X2{Q)]Yi m {Q.)dQ,dx 1 dx2, 

ipjp(lm) actually is the wavelet ipj p in the (/m)-representation. 
Our goal below is to show that $j p (lm) is localized in the (7m)- 
space. 



Using the Fourier representation of S D function, we have 
= J dkidk 2 (B16) 

J ' e ,[i l((1 ) tl+ x 2 ( fl )M yim(!])(|S! 

dk 1 dk24> jp (k 1 ,k2) j dQe iV2Rkcos6 'Y lm (n). 



= J dkidk2i>j p (ki,k 2 ) J ( 



where ipji(ki, k2) is the Fourier transform of the wavelet. In the 
last step of eq.(B16), we used Eq.(B4), ki = k cos 4>, — k sin <j> 
(or k = (fc 2 + kl) 1 ' 2 ) and cos(9' = y/l - |sin0|. Wavelet 
ipji{k\,k2) is localized in fc-space. For a given j = (j'1,.7'2), 
i>ji{ki,k2) is localized in the range 

2n2^ +1 ' 2 /R and 2n2^- 1 ' 2 /R < k 2 < 2n2 j2+1 ^ 2 / R. There- 
fore, for the case ji = j2 = j, the wavelet ipji(ki, fe) is localized 
in the fc-band of 2n2 3 /R < k < 2n2 j+1 /R. 

After the integral on azimuthal angle <j>, eq.(B16) gives 



$jp(lm) = Som^^r^- I dkidk2ipj p (ki,k2) 



2tt 3 / 2 



(B17) 



n I dcosOe 



V2flfct 



P(cos<9)- 



Using the expansion of e 



\/2Sfecos «' 



with P ; - (cos #'), we have 



i>j p (lr 



So, 



]T(2/' + 



v / 2T+7 

" 2tt3/2 

c 

dkidk2tpj p (ki, k2)ji(V2Rk) 
P ( -(cos 6')Pi{cos6)dcos6, 



(B18) 



where ji (x) is spherical Bessel function. First, we take approxima- 
tion cos 9' ~ cos 0, eq.(B 18) gives 



1 n \ x a/2? + 1 .i 
ipj p (lm) = dom 375—* 



dkidk2tpjp{ki,k 2 )ji(V2Rk)(Bl9) 



where ji (x) is the spherical Bessel function. It is known that the 
function ji(x) <C 1, when x < I, it approaches its peak at x ~ 
In, and oscillating around zero when x > In. On the other hand, 
tpjp(ki, fe) is slowly varying with k in the range 2 J / R < k/2n < 
2-' +1 / R. Therefore, the integral on k in eq.(B19) is of non-zero 
mainly in the range of I ~ \/2Rk/n. Thus, from the fc-band of 
2nT jR < k < 2n2 J+1 /R, we have that for large j, 4> jp (lm) is 
of non-zero mainly is in the /-band given by 

2 i+3/2 < j < 2 j+5/2_ ( g 20) 

This result actually is the well known property of the localization of 
wavelet in scale space, regardless that the scale-space is described 
by k- or /-representations. 

Now, let us consider the effect of cos 8' 7^ cos 0. One can 
estimate this effect by the expansion P t r (cos 9') = P ; /(cos6> + 
77) = P;- (cos 0) + ][](rf '/n!)P/, n) (cos60, where 77 = cos6»' - 
cos 6 = VI — sin 6 — cos 6. The derivative of Legendre func- 
tion, P/ n) (cos6l), can be expressed by a summation of P;_ n ... 
P; +n . Therefore, the n 4 ' 1 correction may enlarge the /-band from 
2 i+3/2 -2 J+5/2 to(2 J+3/2 -n)-(2 J+5/2 +n).Since77 < 0.3, we 
need to consider only a few low order corrections, say n < 3. On 
the other hand, for high j (small scales) case, we have 2 J >• 1, or 
2 J » n, and therefore, the r\ correction on the /-band of tpj p (lm) 
is small when j is large. 

One can then conclude that, for different j, the wavelet 
tpj p (lm) consists of a,j m in different /-band, or the /-bands for j 
and j' (j / j') basically are not overlapped. From eq.(B14), the 
coefficients a im in different /-bands are uncorrelated. Thus, the co- 
variance (ejp£j'p>) should be diagonal or quasi-diagonal as 



(tjptj'p') — 0, 



(B21) 



